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We introduce on/off intermittency into a mean field dynamo model by imposing stochastic fluctuations in either the alpha 
effect or through the inclusion of a fluctuating electromotive force. Sufficiently strong small scale fluctuations with time 
scales of the order of 0.3-3 years can produce long term variations in the system on time scales of the order of hundreds of 
years. However, global suppression of magnetic activity in both hemispheres at once was not observed. The variation of the 
magnetic field does not resemble that of the sunspot number, but is more reminiscent of the ^"Be record. The interpretation 
of our results focuses attention on the connection between the level of magnetic activity and the sunspot number, an issue 
that must be elucidated if long term solar effects are to be well understood. 
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1 Introduction 

The sun shows variability on a broad range of time scales, 
from milliseconds to millennia and on to even longer scales. 
Here we are interested in the time scales associated with the 
solar cycle and its grand minima. A salient manifestation 
of the cyclic behavior is seen in the sunspot number, whose 
annual mean oscillates on a scale of eleven years (or twenty- 
two years, if one goes by magnetic polarity variations). The 
term "cycle" is used to describe this oscillation in the sense 
that the sunspot number qualitatively performs the same 
kind of oscillation approximately every eleven years, with 
an amplitude that varies in a way that is reminiscent of some 
chaotic oscillators. However, it has been found that we do 
not have sufficient data to decide whether the global solar 
magnetic variation is chaotic, in the sense that its temporal 
behavior may be that of a low-order deterministic dynami- 
cal system (Spiegel & Wolf 1987). Nevertheless, a chaotic 
oscillator offers a natural way to model various irregularities 
of the solar cycle (Spiegel 1977; Tavakol 1978; Ruzmaikin 
1981). 

Like certain simple chaotic systems, the sun repeats it- 
self magnetically over and over again, but never quite the 
same way twice, much in the manner of simple chaotic os- 
cillators. On the other hand, the simplest chaotic dynamos 
(Allan 1962; Robbins 1979) do not exhibit the kind of 
strong intermittency such as was seen in the Maunder Min- 
imum (Eddy 1978) when the amplitude of the oscillation 
in sunspot numbers went nearly to zero for seventy-five 
years in the time of Newton. This behavior is suggestive 
of the possibility that a strongly intermittent chaotic oscil- 
lator needs to be invoked in trying to understand the so- 
lar oscillation. Oscillators that behave this way are known 
(Spiegel 1981; Fautrelle & Childress 1982) but, when they 
are in the active phase, their variations do not normally re- 
semble those of the sunspot number. 
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In modeling the solar magnetic variability with simple 
oscillators, it is not obvious how to connect the output of 
the models with the sunspot number It may therefore not 
be damning if the variations produced by a model do not re- 
produce even qualitatively the variations seen in the sunspot 
number Nevertheless, in terms of nonlinear lumped mod- 
els, whose behavior is purely temporal, it has been pos- 
sible to produce variations of the cyclic behavior that re- 
semble those seen in the grand minima qualitatively either 
by modulation (Weiss et al. 1984) or intermittency (Piatt 
et al. 1993b). But the solar variability is manifestly spatio- 
temporal and lumped models can at best provide clues to 
the actual processes involved in the solar cycle. Here we de- 
scribe an attempt to go beyond models with purely tempo- 
ral variations to models showing spatio-temporal variations 
that also produce analogues of the grand solar minima. 

An earlier version of this paper was composed in the 
early nineties and some added remarks were inspired by 
discussions at the Enrico Fermi School in Varenna (Cini 
Castagnoli & Provenzale 1997). In the meantime a lot of 
new work has emerged, but we feel that the ideas presented 
here are still relevant. Particularly important has been the 
work of Beer et al. (1998) in demonstrating the persistence 
of a solar activity cycle throughout the time of the Maunder 
Minimum. Such behavior emerged from a purely temporal 
model (Pasquero 1996) as well as from a one-dimensional 
dynamo model (Tobias 1996) with small turbulent magnetic 
Prandtl number (so that the viscous diffusion timescale is 
much longer than that for magnetic diffusion time). Similar 
results — also for small turbulent magnetic Prandtl num- 
bers — have been obtained for two-dimensional models by 
incorporating quenching in the A-effect that drives the dif- 
ferential rotation (Kiiker et al. 1999). Intermittent behavior 
in dynamo models has been studied further by Tworkowski 
et al. (1998) using time-dependent alpha-quenching, and by 
Charbonneau et al. (2005) using algebraic alpha-quenching. 
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However, intermittency can equally well occur in stochas- 
tically forced models (e.g. John et al. 2002). In the fol- 
lowing, we discuss the spatio-temporal variability of sim- 
ilarly forced dynamo models in two dimensions assuming 
axisynmietry. 

2 Spatio-temporal variability 

The spatio-temporal dynamics of the solar cycle as seen in 
spacetime diagrams like the Maunder butterfly diagram sug- 
gest that sohtary waves may play an active part in the solar 
activity cycle (Proctor & Spiegel 1 99 1 ) . A nonUnear version 
of Parker's (1955) dynamo waves (Worledge et al. 1997) or 
solitary waves arising from another overstability could be 
the mechanism of the drifting of the center of activity in 
latitude through the cycle. The finite width of the activity 
zone at any given time suggests that the waves in question 
are themselves confined or guided by a layer of some cor- 
responding thickness. The picture that we adopt here is the 
now conventional one that the layer is the tachocline, though 
such a layer has been variously considered to lie deep in 
the convection zone (DeLuca 1986), just below it (Spiegel 
& Weiss 1980) or occupy the full convective zone (Bran- 
denburg 2005). It might operate on the standard ingredients 
of a dynamo — differential rotation and cyclonic convec- 
tion — from which we here make a model. Since we first 
wrote those lines, the layer has been well studied both ob- 
servationally and theoretically (Hughes et al. 2006). It me- 
diates the transition between the outer differentially rotat- 
ing layers of the convection zone and the inner core with 
its nearly constant angular velocity. It has been renamed the 
tachocline (Spiegel & Zahn 1992) in keeping with its gain 
in respectability. 

Here we use mean field theory (Moffatt 1978; Krause & 
Radler 1980) with the effects of small scale motions sub- 
sumed into turbulent diffusivity and the a-effect. Though 
mean field dynamos may not tell the whole story of the 
solar magnetic fluctuations (e.g. Hoyng 1987; 1988), they 
will suit our purpose of modeUng the intermittency signaled 
by the Maunder Minimum. In modeling solar intermittency, 
we need to be aware that there are several forms of inter- 
mittency that have been isolated in dynamical systems the- 
ory. However, there is a particular one that models the solar 
grand minima quite well (Pasquero 1996) and that has come 
to be called on/off intermittency (Piatt et al. 1993a; see also 
Spiegel 1994). [In fact, this form of intermittency was fash- 
ioned (Spiegel 1981) with grand minima in mind.] 

On/off intermittency is Uke the intermittency detected in 
the output of a probe in a turbulent fluid registering abrupt 
changes as laminar and turbulent fluid regions flow past it. 
One may think of this as a series of bursts or chaotic relax- 
ation oscillations. What characterizes models that have been 
made of this process (e.g. Spiegel 1981; Ott & Chen 1990; 
Pikovsky & Grassberger 1991) is that the (potentially) un- 
stable oscillator performing the cycle is driven to instability 
through couphng to an aperiodic driver that is continuously 



chaotic or stochastic. The role of the driver is to move the 
system into and out of the unstable state. Observations of the 
bursts cannot readily distinguish the chaotic driver from a 
stochastic alternative which works equally well (von Hard- 
enberg et al. 1997), and the number of degrees of freedom 
involved in such an object in the solar case will be hard to 
determine (but see Heagy et al. 1994). Fortunately, the qual- 
itative features of the process do not depend sensitively on 
this difference, and we shall use a stochastic driver here. In 
our present considerations, the action of the driver is meant 
to model the influence of the convective solar dynamo and 
the on/off oscillator represents the tachocline. Here we in- 
troduce on/off intermittency into a working model of the so- 
lar cycle (Rtidiger & Brandenburg 1995) that has not previ- 
ously produced grand minima. Modulational grand minima 
in dynamo models have been found by Tobias. In the case of 
a distributed dynamo (Brandenburg 2005), on/off behavior 
may similarly be produced. 

The effects of stochastic noise on mean field ail dy- 
namos have been studied previously (Choudhuri 1992; 
Moss et al. 1992; Hoyng et al. 1994), mainly to model irreg- 
ularities of the solar cycle on time scales comparable to the 
solar cycle itself or shorter. Moss et al. (1992) suggested that 
a modulation of the solar cycle on longer timescales of the 
order of centuries should also be possible. More recently, 
work along those lines (Schmitt et al. 1996) has introduced 
the on/off intermittency mechanism into a mean-field style 
of dynamo, as had been used already in a lumped model 
of the solar cycle (Piatt et al. 1993b). The present paper is 
similarly based on a relatively reaUstic model of the solar 
dynamo in being spatio-temporal and into which we intro- 
duce the on/off intermittency mechanism. 

The main difference between our work and that of 
Schmitt et al. is that we consider a standard a-effect while 
they introduced a lower cutoff excluding field generation 
below a field strength of 1 kG at the base of the convection 
zone. This kind of filtering is analogous to Dumey's (1995) 
introduction of a critical field value in the tachocline above 
which a magnetic tube is ejected. This also produces grand 
minima in spatio-temporal models. 

Another difference with the calculation of Schmitt et al. 
is that we consider the full sun whereas they studied only 
a single hemisphere and produced one-winged butterflies. 
This relates to another aspect of the problem that our model 
is intended to bring out. It is very difficult to have both hemi- 
spheres go completely inactive for appreciable times. The 
probability of a complete tumoff is exceedingly low and it 
seems unlikely that any propagative intermittency mecha- 
nism with retardation effects could turn off the magnetic 
activity globally in models with large spatial extent. The 
model described here does however produce lowered solar 
activity over large portions of its computational domain and 
therefore gives the kind of reduction in total overall activity 
that is consistent with what was seen in the Maunder Mini- 
mum. 
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The hemispheric asymmetry has been modeled by 
Knobloch et al. (1998; see also Weiss 1993) and this work 
is related to the problem of making an extended system 
demonstrate on/off intermittency. The problem has some- 
thing in common with that involved in laying a rug. Once 
the rug is nailed down, if there is a bump somewhere, there 
seems to be no way to squeeze it out of existence. We have 
observed the same behavior with the dynamo model. If we 
produce a grand minimum locally somewhere in space, we 
find invariably that there is usually some excitation else- 
where. While we can get a whole hemisphere to be quiet 
at one time, the other one may still show activity. The out- 
put of the dynamo model is therefore somewhat like the sun 
in intermission — on the global scale there is usually some 
weak activity somewhere and it does have cycles. [Pasquero 
(1996) has shown how this may be achieved in the purely 
temporal oscillators as well.] As we have learned at the 
Fermi School on the Interaction of the Solar Cycle with 
Terrestrial Activity in June of 1996, certain terrestrial data 
have much more in common with the output of the dynamo 
models than the sunspot number. Moreover, some of these 
terrestrial data are very convincing proxy data for the solar 
activity (Beer et al. 1990; 1994; 1996; Solanki et al. 2004), 
while others may mimic the sunspot number. [There are also 
auroral indicators of solar activity recorded in classical an- 
tiquity (Stothers 1979, Solow 2005).] To a great extent, the 
issue is at heart one of knowing how to compare the output 
of dynamo models to the various measures of solar activity. 

We describe next some particulars of the model itself, 
including the manner in which the fluctuations are intro- 
duced. Readers who are not interested in such detailed infor- 
mation should skip directly to Sect. 4 where we outline the 
main results. In these, we focus on models with a rather high 
noise level exceeding the electromotive force from the a- 
effect by an order of magnitude since we expect the global 
convective dynamo action to be more vigorous than that of 
the tachocline. 

3 The model 

We use a mean-field model of the dynamo action in 
the tachocUne (Riidiger & Brandenburg 1995) with an 
anisotropic a-effect and a turbulent magnetic diffusivity. 
Magnetic buoyancy is also included as it is in the more 
elaborate model of Jiang et al. (2007; see further references 
therein). The prescribed angular velocity is taken from 
the observational findings of heUoseismology (Christensen- 
Dalsgaard & Schou 1988). To obtain a 22 yr magnetic cycle 
period, we introduce a scaling factor in the magnetic diffu- 
sivity of 0.5 and, to get a butterfly diagram with sufficient 
activity at low latitudes, we introduce a suitable latitude de- 
pendence in a. 

In the original model calculations, the value of a was 
typically set at a value approximately twenty times the crit- 
ical value for instability. To produce on/off intermittency 
(Piatt et al. 1993a), in the highly supercritical case, we 



would need very large fluctuations at that value. In this ex- 
ploratory study, we prefer to operate at more modest pa- 
rameter values to avoid the need for large fluctuations in the 
driver. Therefore we choose a slightly subcritical value of 
a and introduce only modest fluctuations in its magnitude 
so that the driver can move it into and out of the unstable 
state easily. For this purpose we introduce a scaling factor 
Cq, in front of certain components of the a-tensor (those 
components which result from the interaction of rotation 
and stratification). For details see the original discussion of 
the model (Riidiger & Brandenburg 1995) where, with = 
1, the resulting toroidal magnetic field is a few kgauss and 
the poloidal field at the surface about 10 gauss. In the cases 
studied here, however, where the dynamo is just marginally 
excited, the generated magnetic field is weak, and could not 
explain the field strength observed in sunspots. This feature 
of the model can be corrected, as the work of Schmitt et al. 
(1996) shows. 

We consider separately two kinds of fluctuation in the 
model: fluctuations in a or in an imposed electromotive 
force S representing the effect of the fluctuations in the 
main solar dynamo. The latter is a more elementary pro- 
cess that does not rely much on dynamo theory. The picture 
here is that in the bulk of the convection zone a small scale 
dynamo operates (e.g. Meneguzzi & Pouquet 1989; Nord- 
lund et al. 1992) producing a magnetic field that is highly 
variable in space and time as represented by £. In the bulk 
of the convection zone the fluctuations are immense, hence 
even spatial and temporal averages (£ ~ {u' x B')) re- 
main fluctuating, albeit on longer scales (Brandenburg et 
al. 2008). The angular brackets refer to ensemble averages 
in principle but, in practice, they are approximated by spa- 
tial and temporal coarse graining averages. The error result- 
ing from this approximation is sometimes interpreted as a 
source of stochastic noise (Hoyng 1987; 1988; 1993; Moss 
et al. 1992; Brandenburg et al. 2008). As we have men- 
tioned, there is ample reason for expecting fluctuations to 
appear in any realistic model. 

We adopt white noise with vanishing mean value and a 
root mean square value of unity. The temporal power spec- 
trum of the noise is flat for frequencies smaller than 2-k/t. 
The value of r determines the time span over which long 
term variability of the resulting mean magnetic field is pos- 
sible. The amplitude of the noise, ng, is measured in terms 
of the rms velocity of the turbulent motions and the lo- 
cal equipartition field strength. Whichever of the fluctuation 
mechanisms apphes — in the emf or in the a effect — the 
procedure is similar. For the a case, for example, we add a 
fraction of the noisy component to the original a. 

4 Results 

4.1 Fluctuating a-effect 

To regulate the strength of the effect of noise, we introduce 
a factor Ca in front of the familiar a term, as described in 
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Fig. 1 Butterfly diagrams of the toroidal field at the base 
of the convection zone for Cq, = 0.02, = 5 and r = 3 yr. 



Sect. 3. With everything else fixed, instability occurs when 
Cq exceeds the critical value Cq"^''-* « 0.03. We begin by 
adopting the subcritical value = 0.02 and we adjust the 
fluctuations in a to have a time scale r = 3 yr. 

If the noise level, , is too low, the magnetic field de- 
cays to zero, as it does when = 0.02. Short noisy bursts in 
a are insufficient to bring the magnetic field to appreciable 
strength. On the other hand, if ria is rather larger than this, 
the magnetic field is almost entirely dominated by the noise 
and cyclic behavior does not occur, except for irregular re- 
versals on the timescale of centuries. Such a case is depicted 
in Fig. 1 , where we show a color-coded representation of the 
toroidal magnetic field at the bottom of the convection zone 
as a function of time and latitude. In the following we refer 
to such representations as butterfly diagrams. A very simi- 
lar result, with reversals on a long time scale, is found even 
when the non-random component of a is absent altogether 
provided there is global shear in the tachocline, as there is 
in accretion disks (Vishniac & Brandenburg 1997). 

To get sufficiently large magnetic field strengths, we fo- 
cus attention on models with = 1. With r = 3yr we 
find long term variability on a time scale of 200-500 yr; see 
Fig. 2. We associate this variability with grand minima (and 
maxima) as discussed below in Sect. 5 although as noted al- 
ready, there is never a complete turning off of the cycle at 
all locations at once. 

Even if r is decreased to a value of 0.3 yr, say, long 
term variability still occurs, but individual cycles may vary 
significantly in amplitude as seen in Fig. 3. Such behavior 
is found only if the fluctuations of a are sufficiently larger 
than the average value; here a factor of five is needed. Indi- 
vidual fluctuations of a can still be much larger because we 
have assumed an exponential distribution for the probability 
density of a fluctuations. 



4.2 Magnetic noise 

If we introduce an external stochastic emf, as described in 
Sect. 3, we find a qualitatively similar behavior to that with 
fluctuating a; see Fig. 4. Fluctuations of various kinds can 
evidently produce long time intermittency. 

The model shows two distinct activity waves, one mi- 
grating equatorward and the other poleward. The two waves 
seem to be modulated independently, in each hemisphere. 
So do the modulations in the two hemispheres seem to be 
only weakly coupled, with little tendency to form a dipole 
structure. It may be that the approximate antisymmetry of 
the toroidal magnetic field of the sun, suggested by Hale's 
polarity law, may not be a very stable feature, and that 
other types of (a)symmetry might have occurred in the past. 
Other, more regular parity variations of the magnetic field 
have previously been seen in nonlinear models (Branden- 
burg et al. 1989a,b; 1990; Jennings & Weiss 1991; Sokoloff 
& Nesme-Ribes 1994). 

5 Interpretation 

In the model described here, the seat of the solar activity cy- 
cle is in the solar tachocline, the layer that matches the dif- 
ferential rotation of the solar convection zone to the (nearly) 
rigid rotation of the inner sun (Hughes et al. 2006). While 
the precise hydromagnetic process that drives the activity 
has not yet been securely identified, we have assumed that 
the process may be modeled as an acu dynamo for the pur- 
pose of exploring the cause of the grand minima of solar 
activity. However, as we have already implied, any of sev- 
eral overstabilities might equally serve our purposes. 

As we have tacitly assumed, the differential rotation in 
the tachocline is likely to be significant in such instabilities. 
In particular it is likely to give rise to a toroidal field. Given 
a suitable depth dependence of the strength of this field, 
an instability driven by magnetic buoyancy (Parker 1979) 
may arise and, especially in the presence of a stable density 
stratification, it may drive waves of excitation (Proctor & 
Spiegel 1991). Another possible driver for such waves may 
be magnetorotational instability (Balbus & Hawley 1998). 
Parfrey & Menou (2007) have studied the local instabil- 
ity of the tachocline and found instability at high latitudes 
and stability at low latitudes. However, their calculation did 
not include a toroidal field and, as Knob loch (1992) has 
observed, an azimuthal field can cause a Hopf bifurcation 
in MRJ. As we may reasonably expect to find a toroidal 
field in the tachocline, this overstabiUty is another mecha- 
nism that could perhaps engender waves whose description 
would typically be by the complex Ginzburg-Landau equa- 
tion (Aranson & Kramer 2002). A related calculation is de- 
scribed by Kitchatinov and Riidiger (2007); see also Cally 
(2003), Oilman et al. (2007), Rudiger & Kitchatinov (2007), 
and Zahn et al. (2007). 

Given that the tachocline is poised to drive hydromag- 
netic activity, by dynamo processes or otherwise, we have 
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Fig. 2 Butterfly diagrams of the toroidal field at the base of the convection zone for Ca = 0.02, Ua — 5 and r = Syr. 
Butterfly diagrams of the toroidal field at the base of the convection zone for = 1, ria = 0.2 and t = 3 yr. 




Fig. 3 Butterfly diagrams of the toroidal field at the base 
of the convection zone for Ca — I, Ua = 5 and t = 0.3 yr. 




Fig. 4 Butterfly diagrams of the toroidal field at the base 
of the convection zone for Cq = 1, = 1 and r = 3 yr. 



also assumed that the ambience is highly fluctuating. Here 
we have in mind that the main dynamo in the sun is a global 



convective dynamo. It is this process that is presumably re- 
sponsible for the magnetic carpet (Simon et al. 2001) in a 
continuous process that is thought to produce rapidly fluc- 
tuating fields of moderate strength (Priest et al. 2002). In 
view of the complications in the theory of this process, we 
have simply modeled the highly fluctuating fields in of the 
solar convection zone as a stochastic process. In the mech- 
anism of on/off intermittency that we have introduced here, 
the fluctuations produced by the main dynamo move the 
tachocline into and out of states of instantaneous oversta- 
bility. This mechanism produces, as we have noted, longer 
response times in the model tachocline than the time scales 
of the fluctuations that give rise to more intense field con- 
centrations. [We have seen a similar kind of symbiosis in 
the simulations of the geodynamo by Glatzmaier & Roberts 
(1995)]. The details of the process are complicated but, in 
gross, when the local field fluctuations produced are too 
weak, the grouping of flux ropes into a sunspot field is not 
achieved, even though much of the normal activity contin- 
ues in the convection zone. To illustrate how such details 
may relate to the observed spot numbers, we have mapped 
the activity to the butterfly diagram with a particular func- 
tional of the field in the tachocline for the present purposes. 

Like the sun, the model we have studied here mani- 
fests spatio-temporal intermittency but the magnetic activ- 
ity almost never turns off everywhere. That is, our results 
suggest that simple dynamo models, with either fluctuating 
emfs and/or with magnetic noise injected into the bulk of the 
convection zone, can produce intermittency on sufficiently 
long time scales, but they do not switch off globally over the 
whole sun, in both hemispheres at once. It seems likely that 
this feature is inherent in all models with propagative behav- 
ior and that it cannot be expected that the cycle-producing 
processes of the sun switches off everywhere for several cy- 
cles as might be imagined on the basis of certain lumped 
models (such as Piatt et al. 1993b). When we first produced 
the present results, we thought this feature was a deficiency. 
But as Ribes 8c Nesme-Ribes (1993) have reported, even 



5 



A. Brandenburg & E. A. Spiegel: Modeling a Maunder Minimum 
Astron. Nachr. /AN (0000) 




500 550 600 650 



700 750 800 



to 

(0 




Fig. 5 Time series of | _B | (dotted line) and the activity pa- 
rameter R (solid line) for the same run as in Fig. 2. Note that 
the scale for increases downwards, in order to mimic the 
approximate anticorrelation between the ^°Be data and the 
sunspot number. 

during the Maunder Minimum, a weak solar cycle contin- 
ued. Thus a globally depleted activity level is more like what 
is wanted and we have gotten that from the model in keep- 
ing with the historical records as interpreted by Ribes and 
Nesme-Ribes. 

There still remains the need for an additional physical 
feature that relates to the production of sunspots, or more 
specifically, strong flux tubes of sufficiently large cross- 
section. That is the message of the procedures of Durney 
(1995) and Schmitt et al. (1996). A simple way of formulat- 
ing this problem is to say that the sunspot number, which is 
a global parameter, is, as already mentioned, a functional of 
the various fields produced in the model. To get that func- 
tional, we need to operate with an expUcit sunspot produc- 
tion mechanism. 

Let R{t) be the average magnetic field strength between 
±10° and ±30° latitude, but allow only those areas where 
the field exceeds the root mean square value by 30 percent to 
contribute to this average. This quantity is plotted in Fig. 5. 
Note that R{t) shows periods with almost vanishing mag- 
netic activity as in a Maunder Minimum. The grand min- 
ima on this interpretation of the output of the models result 
from what may be regarded as lean cycles magnetically. If 
there are magnetic droughts, they are only local, not global, 
though the magnetic means may be low. In that sense, the 
sunspot number, though valuable for having focused our at- 
tention on an interesting feature of the magnetodynamics, 
may in some ways be a misleading indicator of what is hap- 
pening overall. 

This view is supported by studies of other indicators of 
solar activity than the sunspot number. The most complete 
records of variations that may result from solar activity fluc- 
tuations are found in the ^"Be records from ice cores (Beer 
et al. 1990; 1994; 1996). The i°Be variations are plausibly 
attributed to modulation of the cosmic ray flux by the mag- 



Fig. 6 "Be data (solid Une) together with the sunspot 
number (dotted line), as provided by Dr. Jiirg Beer. The data 
are from a shallow core (300 m) drilled at Dye 3, Greenland, 
in 1986. The data were filtered using a spectral filter with a 
cut-off of 6 years and interpolated using a cubic spline. The 
younger part (1783-1985) is pubhshed in (Beer et al. 1990), 
and the whole record appeared in Beer et al. (1994). 

netic field in the solar wind (Beer et al. 1996). In Fig. 6 
we reproduce data kindly provided by Dr. Jiirg Beer com- 
paring the ^°Be records with the sunspot number for nearly 
four centuries. We see that a cyclic modulation is very much 
in evidence during the Maunder Minimum, with a rather 
modest reduction in its amplitude compared to that of the 
sunspot number. 

The situation as brought out in the cited papers of Beer 
et al. is that various measures of solar activity are not per- 
fectly correlated. If we were to think of the output of a 
model solar dynamo as we would one of these other mea- 
sures, we would not be surprised that it does not necessarily 
represent all of them faithfully. Unfortunately, there is as yet 
no clear theoretical indication which one of them a model 
should most closely represent. If, in a grand minimum, there 
is high magnetic activity in only one hemisphere, it is not 
unreasonable that the modulations of ^°Be should continue 
with reasonable strength. The question for the theory then 
is what is the relationship between the variations produced 
by the models and (say) the sunspot number. This is partic- 
ularly difficult since the sunspots seem to be produced well 
within the sun and not at the surface, whereas some of the 
other activity measures are no doubt superficiaUy produced. 
Even if we do have a good model of the cyclic mechanism, 
we also need to understand how sunspots form and surface 
before we can predict their number. 

Nor is it unimportant to try to predict the sunspot num- 
ber, or something akin to it, for it is this quantity that seems 
to be connected to some climatological variations. The most 
striking evidence of this is the discovery in tree ring data 
(Douglass 1927) that "the sunspot curve flattens out in a 
striking manner ... from 1670 or 1680 to 1727." This dis- 
covery was made by A.E. Douglas before he had "received 
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a letter from Professor E. Maunder ... calling attention to 

the prolonged dearth of sunspots between 1645 and 1715." 
As Douglass and others have argued, variations in tree ring 
thickness in turn are connected to rainfall, so it is not an idle 
project to try to understand how to go from the workings of 
a solar dynamo to the manufacture of sunspots. 

In summary, solar activity waves at high and low lati- 
tudes and in the two hemispheres, lead somewhat indepen- 
dent, weakly correlated hves and fluctuate separately under 
the influence of noise. Spatial variations of the solar cycle 
during the Maunder Minimum are not an immediate indi- 
cator of the sunspot number and, if the sunspots have their 
origin well beneath the solar surface, we must go another 
step in the discussion before we obtain results that are fully 
consistent with historic records of sunspots as summarized 
by Ribes and Nesme-Ribes (1993). It is not at all obvious 
what we may conclude about subconvective activity from 
observed surface activity. 
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